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ABSTRACT 

+ 
The  capture  cross  section  for  the  formation  of  H   in  the 

+  +  . 

reaction  H   +  CH,  -•  H„  +  CH0  was  measured  at  incident  proton 
-4—23 

energies  of  JO,    85,  100,  150  and  200  eV  and  covering  the  scat- 
tering angles  of  43   to  49-5   (lab  coordinates).   At  100  eV  and 

below  the  curve  of  the  cross  section  versus  angle  shows  a  sharp 

;  /■  o  .   .  ... 

peak  at  about  4°   whose  position  approaches  the  theoretical  limit 

of  46.9   with  increasing  energy.   Above  100  eV  the  peak  was  too 

small  to  be  observed  and  only  an  upper  limit  can  be  placed  on 

the  value  of  the  cross  section.   Typical  values  of  the  total 

-21    2  -22 

cross  section  are  2.0  x  10     cm   at  70  eV  and  7.6  x  10   'at 

100  eV .   The  magnitude  and  energy  dependence  of  the  cross  section 

as  well  as  the  angular  position  of  the  peak  all  are  in  essential 

agreement  with  the  classical  theory  of  ion-molecule  rearrangement 

collisions  proposed  by  Bates,  Cook  and  Smith. 
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I.   INTRODUCTION 

The  correct  description  of  complex  atomic  and  molecular  inter- 
actions is  one  of  the  major  achievements  of  modern  quantum  theory. 
While  there  appears  to  be  little  doubt  about  the  correctness  of 
the  theory  in  principle,  it  has  become  evident  in  practice  that 
the  mathematical  complexity  is  such  that  only  the  very  simplest 
cases  can  be  treated  with  anything  approaching  complete  rigor. 
Therefore  many  approximations  are  employed  to  simplify  these 
mathematical  computations.   The  question  then,  of  the  accuracy, 
nature  and  range  of  applicability  of  these  approximations  is  of 
considerable  practical  importance.   In  general  these  questions  can 
only  be  answered  by  comparison  of  the  experimental  results  and  the 
pr edict ions  of  a  specific  approximate  calculation. 

This  paper  presents  the  results  of  a  series  of  measurements 
which  confirm  the  predictions  of  a  classical  theory  of  ion- 
molecule  rearrangement  collisions  at  high  impact  energies  as  pro- 
posed by  Bates,  Cook  and  Smith  [  1]  .   This  theory  uses  an  impulse 
type  approximation  to  obtain  the  cross  section  for  the  capture  of 
a  light  atom  or  ion  from  a  target  molecule  by  a  fast  moving  pro- 
jectile.  Applying  this  theory  to  the  reaction: 

H+  +  CH,  -  H*  +  CH  ( 1 ) 

,  ,       o 
the   cross    section    is   predicted   to   show   a    sharp   peak  at    l+o .  9      with 

an   upper    limit    on    the   magnitude    of    the   cross    section    of 

i  -2  0        2 

I.q.  x    10  cm      at    100    eV    incident   proton    energy.       It    is    further 
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predicted  that  the  magnitude  of  the  cross  section  should  decrease 
rapidly  with  increasing  proton  energy,  approaching  an  energy  de- 
pendence of  E  "    asymptomat ically  at  high  energies  (above  500  eV) . 
The  range  of  validity  of  the  theory  extends  from  about  50  to  800  eV 
when  applied  to  reaction  (1). 

Our  measurements  were  carried  out  at  JO,    85,  100,  150,  and 
200  eV,  covering  the  angles  of  scatter  from  43   to  49.5  .   We 
found  essential  agreement  between  the  predictions  of  the  theory 
and  the  experimental  data,  indicating  that  the  classical  impulse 
approximation  proposed  by  Bates,  Cook  and  Smith  is  a  valid  model 
for  reaction  (1)  in  the  energy  region  considered. 
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II.   THEORY 

A.   CLASSICAL  THEORY  OF  ION-MOLECULE  REARRANGEMENT  COLLISIONS 

A  classical  theory  of  ion-molecule  rearrangement  collisions 
at  high  impact  energies  has  been  proposed  by  Bates,  Cook  and 
Smith  [ l] .   The  impact  energies  are  assumed  to  be  high  enough 
so  polarization  forces  and  chemical  binding  energies  of  the 
colliding  molecules  can  be  ignored.   The  theory  applies  to  various 
ion  molecule  rearrangement  reactions  of  the  type: 

X'X.  +YZ,  Z'  -  X1  +  X.Y  +  (Z,  +Z')  (1) 
1      k      —      1     v  k      '  v  ' 

where  X.  and  Y  are  simple  atoms  or  ions.  The  bar  indicates  which 
ions  are  fast  in  the  laboratory  coordinate  system.  An  example  of 
a  rearrangement  reaction  of  this  type  is: 

H+  +  CH,  -  H^  +  CH  (2) 

The  basic  assumption  made  by  Bates,  Cook  and  Smith,  is  that 
process  (1)  may  be  described  by  a  classical  impulse  approximation 
similar  to  that  developed  by  Thomas  [2],  for  the  description  of 

electron  capture,  where  each  composite  system  is  regarded  as  a 

+         + 
loose  cluster  of  atoms  and  ions.   (Eg  H   =  H  +  H   and 

CH,  =  C  +  3H  +  H) . 

According  to  this  model,  we  can  consider  the  reaction  to  occur 

in  the  following  manner:  consider  a  particle  of  mass  M   moving  with 

— ♦ 
velocity  v   through  molecules  of  loosely  bound  particles  of  mass 

M2  and  M     Figure  1  shows  the  sequence  of  collisions  that  lead  to 
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M. 


E 


FIGURE 


THE      CAPTURE      MECHANISM 
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the  capture  of  M   by  M  .   First  mass  M   collides  with  mass  M_ . 

—4 

This  causes  M   to  be  scattered  at  an  angle  9'  and  velocity  v' 

while  M   recoils  at  an  angle  9'  with  velocity  v ' .   Mass  M   then 

suffers  a  second  binary  collision  with  mass  M,  .   This  causes  M„ 

k  2 

— ♦ 

to  scatter  at  an  angle  9"  with  velocity  v". 

Now  if:  (a)  0£  ^  9^  +  Q^\  [3) 

or  vj  =  v" 
And  if:  (b)  v*  <»  v" 

then  masses   M      and  M      may   have  a    relative   energy    of  motion   below 
that   required  for    separation,    and  so,    the   two  particles    may   com- 
bine and  proceed   on   as    one.       If   the  particle  with   mass   M      is    a 
composite  particle   then    its    disruption    can   be   avoided    if 

v1  <«  v      .  ^ 

11 

In  the  particular  case  of: 
M2  =  Ml  «  Mk 

condition  (3)  is  satisfied  if  9'  ^«  45  and  9"  %f  90  (see  Section 
II  of  this  Chapter),  but  then  condition  (4)  is  violated.  In  the 
case  of  the  reaction: 

H+  +  CH,  -  H*  +  CH„ 

-       4—2      3 

where  M      =   M       ,    M_    =   M      and  M      =  M         ,    the   violation    of    condition 
1  H+         AW  k  CH_ 

3 

(4)  is  not  relevant  since  M..  is  not  a  composite  particle. 

Consider  now,  the  first  binary  collision  between  M   and  M„  as 

shown  in  Fig.  2.   The  probability  of  scattering  M?  into  the  solid 

— »      _«     -♦ 
angle  dft(9')with  velocity  between  v'  and  v'  +•  dv'  is; 

q  =  ff12(©J)  dQ(9^)  (5) 


15 


M 


v 


da  (02) 


FIGURE     2        THE  FIRST    BINARY  COLLISION 


FIGURE     3.      THE  SECOND   BINARY    COLLISION 
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where  o       ( 0 ' )     is    the   appropriate    differential    scattering    cross 
section.    Since    the   solid   angle    is    the   cone    of    semi-angle   0'    within 


d9'    then 


dQ(O')  =    2n  sin   0'    d0' 


and   equation    (5)    becomes: 


q    =    2ttct12(Q2)    sin    °2    dQ2     "  (6) 

Consider    now,    the    second  binary    collision    between    M0    and  M      as 
shown    in    Fig.    3-      The  probability    of  M      having   an    impact   parameter 
between   C    and  f +   d f     and   an   azimuthal    angle    between  0    and  0    +    d0 
at    a    distance   r    from  M,     is    given    by: 

P  =  eae-^%  (?) 

4tt  r 

where  r  is  the  M  -M  equilibrium  internuclear  separation.  But 
particles  incident  on  M,  with  impact  parameter  (>  and  azimuthal 
angle  0  are  scattered  into  the  solid  angle  dr(9")  at  9''. 


Hence:         ^d^=a2k(Q2)  dQ(G2) 


where  CT?w(ep)  i-s  the  appropriate  differential  scattering  cross 

section . 

But:  dQ(0")  =  sin  G"  dG" 

Hence   equation    (7)    becomes: 

p    =  ^2k(©2')    sin   0"    d02-    -5&_  (8) 

4rr  r 

— »    — * 
Now,  for  M   to  capture  M? ,  v'  -  v"  the  relative  velocity  of  M 

and  M  must  be  contained  within  a  volume  of  velocity  space 
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determined  by  D,  the  mutual  affinity  of  these  two  systems.   In 

the  high  velocity  limit,  Bates,  Cook  and  Smith  [l],  state  this 

requirement  as: 

3/2 
(v^)2  sin  0^  d0^  dJ/)  dv^'  =  \|tt  (£5)  (9) 

where    \  =  %  and  u,  =  ^M   in  the  case  of  M   =  M  <  M  .      (1°) 
Substituting  equation  (9)  into  equation  (8)  we  get 

»2k(Q2>  ,2*y2     X ,  s 

P  =  2 ("IT )     2 •  (1X) 

3r  (v£)  dv^ 

Now,  the  capture  cross  section  arising  from  the  M  -M„-M, 
sequence  of  binary  collisions  is: 

Q(M]_-M2-Mk)  =  Vpq  (12) 

where  V  is  a  dimens ionless  factor,  less  than,  or  equal  to  unity, 
which  allows  for  the  fact  that  M  and  U  may  approach  each  other 
such  that  in  the  resulting  state,  the  affinity  is  less  than  that 
assumed. 

Substituting  equations  (6)  and  (11)  into  (12)  one  gets: 

2TT  <J       (0<)a       (8«)  /2D\3/2  V2 

Q(M    -M2-Mk)    =  V*    !      2  sir,    9-    d0|/^  -§-.    (13) 

3r  \      2  j  2 

Changing    from  Lab.    to   Center    of  Mass      coordinates,    one    gets    for 

the   case    of   Mn     =  M_   «    M, 
12  k 

l6w  o      (90°)  "2k(90°)      /  2D\3/2 
Q(m1-m2-m1c)   =  v i2_ Si 1^  (14) 


where  cr    ^    and  o"         denote    the   differential    scattering    cross 
12  2k 

sections    in    the   center    of   mass    coordinate   systems.       This 
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differential    cross    section    is   at   an    energy    of   relative  motion    of 

2 

\   M  v   and  the  affinity  D  is  at  an  M„-M   internuclear  distance 

of  {Ti.   To  obtain  the  total  cross  section  for  capture,  Q,  one 
now  sums  over  all  binary  collision  sequences  which  lead  to  the 
same  final  capture  process.   For  example,  the  total  cross  section 
for  process  (2)  is 

Q  =  4Q(H-H-C)  +  4Q(H-H+-C). 

This  represents  the  two  possible  collision  sequences  shown  below: 


H 


Oh  Oh 

The    (H   -   H   -    C)    Sequence 


H 


hO 


hO 


The    (H-   H      -   C)    Sequence 


Solution    of   Rutherford's    expression  for    the  differential 

scattering   cross    section   when   v      is    high   enough   for    the  atomic 

field   to  be  Couloumbic,    shows    that: 

-11 


Q(Ml-M2-Mk)   ~   vx 


(15) 


The  range   of   validity   of   (14)    does    not   extend    indefinitely 
as    v.    increases.      Two  requirements   must   be  satisfied: 
(i)   The   De- Broglie  wavelength   limitations, 
(ii)    The  Heisenberg   uncertainty   relations. 


19 


Bates,  Cook  and  Smith  [l],  computed  the  required  differential 
scattering  cross  sections  a   ( 90  )  and  a   (90  )  using  classical 
methods  and  the  analytic  representation  of  the  relevant  Har ti ee 
potentials  as  given  by  Byatt  [ 3] •   From  the  results  of  these 
computations  listed  in  Ref.  1,  the  cross  section  for  the  re- 
arrangement reaction: 

H+  +  CH,  -  H*  +  CH0 

is  determined  to  be  as  shown  in  Fig.  4. 

B.   THE  H+  +  CH,  -  H*  +  CH„  REACTION  KINETICS 

-       H         —  2      3 

The  rearrangement  mechanism  for  the  reaction  is  shown  in 

+ 
Fig.  1,  where  M   represents  the  H  ,  M„  the  H  and  M   represents 

the  CH  .   To  determine  the  exact  value  of  the  scattering  angle  0' 

one  has  to  consider  the  kinematics  of  the  arrangement  process. 

If  the  speed  v   is  assumed  to  be  large  enough  so  the  binding 
energy  of  M   to  M,  can  be  ignored,  then  the  conservation  laws 
can  be  applied  in  the  usual  manner. 

Consider  the  first  binary  collision  of  the  reaction  shown  in 
Fig.  1.   From  conservation  of  energy  and  momentum  we  have: 

^v2  =  ^(vj)2  +  fc«2(vp2  (16) 

M  v   =  M  v1  cos  0|  +  M  v'  cos  9^  (17) 

0  =  M  v*  sin  0J  -  M  v'  sin  0'  (18) 


Now:      M   =  M   =  M 
1     2.  H 

Hence:  the  above  equations  become 

v2  =    (v|)2  +  (v^)2  (19) 
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FIGURE  4 
PREDICTED  CROSS  SECTION  FOR 
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v   =  v*  cos  0'  +  v'  cos  ©'  (20) 

0   =  vj  sin  0'  -  v«  sin  ©'  (21) 

Now  0'  can  be  eliminated  between  (20)  and  (21)  by  rearranging, 
squaring,  and  adding  the  two  equations  to  yield: 

v2  +  (vj)2-  2v1v|  cos  ©j  =  (v^)2  . 
But  from  (19):     (v^)2  =  v2  -  (vj)2 


Thus 


v2  +  (v|)2-  Zv^^    cos  Oj  =  v2  -  (vj)2 


2 
Hence:  (V-J)   =  v1v1'  cos  e-! 


.'.  vj  =  v   cos  0'  (22) 

Substituting  equation  (22)  back  into  equation  (19)  we  get: 

V1    =   v1  (cos  ©j)   +  (v^) 

or:  v2  (1-cos2  0|)  =  (v^)2 

.*.  v^  =  vlSin  ©J  (23) 

Equations  (22)  and  (23)  give  us  v*  and  v'  in  terms  of  v   and 
0* .   Now  we  wish  to  obtain  v'  and  v'  in  terms  of  v   and  0 ' .   Thus 
we  arrange  equation  (20)  and  (21)  such  that  when  we  square  and 
add  them,  we  eliminate  ©'  and  so  obtain: 

v2  +  (v^)2  -  2v±v^    cos  0£  =  (v{)2  . 
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2  2  2 

But    from    (19)    we   have:  (v*  )       =    v       -    ( v '  ) 

2  2  2  2 

Thus:  v       +    (v£)       -    2viv2    COS    °2    =    Vl    "    (VP 

Hence:  v'    =    v      cos    0'     .  (24) 

Substituting    equation    (24)    back    into    equation    (19)    gives: 

2  /       .  s2  2  2     rvt 

v      =    (v|)       +    v      cos      Qx2 

or:  v^(l-cos2Q,p    =    (v')2 

,\       v|    =    v±    sin   0^     .  (25) 

Consider    now,    the   second   binary   collision.      Again    ignoring    the 
binding    energies    we    can   write    the    conservation    of    energy    and 
momentum    equations    as: 

h  M2(v^)2    =   h  M2(v£)2    *   h  Mk    v2  (26) 


M„v'    =  M_v    "    cos    0"    +   M.  v.     cos    © 
2222  2  k    k  k  ( 2/) 


0    =   M„v"    sin    0"    -    M, v,     sin    0,  (28) 

2    2  2  k    k  k  v        ' 


Now:  M      =    M      =   M 

1  2  H 


and:  Mk    =   M 


Thus    the   above   three    equations    reduce    to: 

M 

<V2>       =    *V2»2    +    ST      Vk  <29> 

H 

Mk 
(v.)        =    vm    cos     Qm     +_    ^    cos    gk  (30) 

H 
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Mk 

0  =  v"  sin  0"  -  —   v,  sin  0,   . 
2       2    M     k       k 
H 


(31 


As  we  are  not  interested  in  0,  we  can  eliminate  it  by 

k 

suitably  rearranging  equations  (30)  and  (31)  and  then  squaring 
and  adding.   Thus  we  get: 


(v^)2  +  (v^)2  -  2v«  v2"  cos  Q2={^)     v" 


(32) 


But    from    equation    (29): 


Mk         2    -    /     n2         /       „x2 

»r  vk  -  (v2)    -  (v2 ")  • 

H 


Thus    (32)    becomes 


(v.p2    +    (v'p2 


M  M 

2v.v,  cos  „„  =  _t  (v,,  .  Js  (v„, 

H  H 


or 


(V2>2    f1"    fijO     +    (V2*)2    (1+   ^    J"    2    V2V2    C°S    Q2    =    ° 


M 


Multiplying    through    by 


H 


M       +   M         v 
k  H    \    2 


—  yields: 


MH-Mk 
MH+Mk 


+      ~         -    2 


M 


H 


VMk 


COS     0" 


=    0 


(33) 


M. 


Now,    let:       a 


H 


MR+MH 


and   b    = 


Mk-MH 
Mk+MH 


Then    (33)    becomes 


-    2 


2/ 


a    cos    0" 


b    =   0 


(34) 
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But  now,  for  capture,  we  require  that;  0"  «s  0'  +  0' 


thus:   cos  0"  ^  cos  (01  +  0 ' ) 


ie. ,    cos    0"  x    cos    0'    cos    0'     -    sin    0'    sin    0'       and    using    (24)    and 


(25)    to    eliminate   0'    we   have: 


fV2\  [Vl 

cos    0"  ^   cos    0'  (   —      -    sin   0'      — 

2  llvJ  l\v1 


Hence   equation    (34)    becomes: 


-    2a 


v"    v' 

r    —       COS     0* 

V         V 

2       1 


v"   v. 


-|  —   sin    0' 

V  f     V 

2      1 


b  =   0 


or:       (v») 


2aV2 


M,,   • 


v"v 


22  n. 
cos     0,' 

v,  1 


v"v 


—   sin    0« 

Vl  X 


b(v^)    =  o  .  (35: 


But    now,    the   second   capture   criteria    is    that:    v 


»  ^    ,t  i 


But,    from   equation    (22)    we    have:    v' 


v      cos    0 ' 


Hence,    we   require    for    capture; 


v   cos  0' 


(36) 


Substituting  equations  (22)  and  (36)  into  (35)  we  get: 

v2  cos2  0'  -  2av'  [v'  cos20»-  v  sin  0'  cos2© ' ] -b( v ' )2=  0 

A.  J.  <s  £*  _L  J.  X  X.  £* 


Using    (23)    to   eliminate    v'    from    this    equation   we   get: 


2         2  2  2  2         2 

v   cos    0*     -    2av    sin0'    [v    sinS'cos    0r-    v   sinO'cos    0']-bv    sin    0' 


=  0 


2     2  2     2 

which  reduces  to;   v   cos   0'  -  b  v   sin   0'  =  0. 


Hence,  for  capture  to  occur  we  must  have: 


tan  6;  =   - 


M   +  M 
k    H 

M   -  M 
k     H 


(37) 
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Now  M.  =  M_„   =  15.03506  and  M„  =  1.00797  amu.   Thus  the  angle  at 

which  capture  occurs  is: 

9|  =  46.926°  (38) 

Hence  we  expect  the  capture  cross  section  to  show  a  pronounced 


peak  at  this  angle.   From  Fig,  4?  the  magnitude  of  the  cro 


ss 


-20    2 
section  is  predicted  to  be  approximately  1.4  x  10     cm   if  the 

energy  of  the  incident  proton  beam  is  100  eV . 


C.   THE  RANGE  OF  VALIDITY  OF  THE  CLASSICAL  DESCRIPTION 

A  lower  limit  on  the  range  of  validity  of  the  classical 
treatment  of  the  rearrangement  collision  is  imposed  by  the 
assumption  made  in  the  theoretical  development  that  the  energies 
involved  in  the  collision  are  much  larger  than  the  binding 
energies.   Hence,  when  the  energy  of  the  incident  proton  beam  E 
is  less  than,  or  equal  to  about  10  times  the  binding  energy  we 
expect  the  peak  in  the  cross  section  to  become  less  pronounced 
and  more  spread  out.   Thus  the  theory  is  valid  when  E   ^  50  eV . 

An  upper  limit  on  the  range  of  validity  of  the  theory  is 
determined  by  the  two  requirements: 

(i)  The  Di  Broghlie  Wavelength  K    of  any  particle  involved 
must  be  much  less  than  the  smallest  distance  involved 
in  the  collision  S. 
Hence  we  require:  X   «    S  (39) 

(ii)  The  uncertainty  in  energy  introduced  by  specifying  the 

transverse  position  of  the  proton  must  be  much  less  than 
the  binding  energy  of  the  final  product.   Now,  from  (39) 
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we  know  the  uncertainty  in  the  proton  position  X  is  less 
than  the  distance  of  closest  approach  S. 
But  we  can  write  the  uncertainty  relation: 

AxAp  =2  (40) 

as  SAp  »  — 

Thus  kp   >^   —  (41) 

2        2 

H     8MUS 

H 

We  then  demand  that  the  binding  energy  D  be  greater  than  AE,  or: 

d  »   -  (43) 

8mhs 

To  find  what  restrictions  equations  (39)  and  (43)  place  on 
the  energy  of  the  incident  proton  beam  we  have  to  relate  the 
distance  of  closest  approach  S  to  the  energy  of  the  proton  beam 
E  .   This  can  be  done  quite  readily  in  the  case  of  coulomb 

collisions.   By  considering  the  collisions  to  be  coulomb  it  will 

C 

be  shown  that:        1 

and  A.  = 


FT 


where  C   and  C      are  constants.   Because  S  decreases  faster  than  \ 
as  E..  increases  we  see  that  equation  (39)  will  indeed  place  an 
upper  limit  on  the  energy  E  . 
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In  a  Coulomb  collision  the  point  of  closest  approach  is 
yielded  by  solving  the  equation: 

i   -  ^  "  ^  =  o  (44) 

s  ui 

(see    equation    3_4~3    °f    Ref.    4) 

2 
U    ,    the   Center    of  Mass    energy,    can    be   expressed  as   U      =   %  u    v 

where  u    is    the   reduced  mass    and   v       ,     is    the   relative    speed   of 

r  el  r 

approach.   0(S),  the  potential  energy  at  the  distance  of  closest 

approach,  is: 

0(s)  .  (Ze)(Z'e') 

The    impact    parameter    b    is    given    by: 

»..     <Ze><z'e'>     cot-f  («) 

U     V 

rel 
(see  equation  3-8-8  of  Ref.  4) 

Where  (R)  is  the  Center  of  Mass  scattering  angle. 
Now  we  have  two  collisions  to  consider: 

(i)  In  the  first  binary  collision  (H  -H)  we  have  (see  Fig.  1) 

M„ 
®  =  90°,  u,  =  -f   ,    vrel  =  v1  (46) 

(ii)    In    the   second   binary    collision    (H-C)    we   have    (see    Fig.    1) 
(g)    ^   qm    =    90°,    a   «  MH    ,    vrel    =    y£  (47) 

and   since      E_   ^  h  E,     we    have   v'    x  .=■  vn     . 
2  1  2^2       1 

Making    these    substitutions    in    (45)    we    get: 

b  =|-  (48) 

H+-    H         El 

Vc     =   12  i^  tW 


i  2 

where:  E      =   *g  M      v    . 

1  Hi 

Hence    for    the   H    -    H   collision,    equation    (  45)    yields 

2 
Which   gives:  S(H  +  -    H)    **   2.4   §-    .  (50) 

El 

Similar ily ,    for    the   H-C    collision   we   get: 

s2    .  lis!      s   -   i%^   =   0 

2 
which  gives:       S(H-C)  «  7-4  |~  •  (51) 

+         .  . 
Hence  the  H  -  H  collision  has  the  smallest  distance  of  closest 

approach  and  hence  equation  (50)  will  establish  the  upper  limits 

of  validity  of  the  classical  approximation. 

Consider  first  the  restrictions  imposed  by  equation  (  39)  •   We 

require  that:      X  «  S 


h      «  2.4  %- 


f«H 


I?  El 


.  E  «  2.32  x  10    Ergs  (52) 

•  •   x 

or     E  «  1.45  x  10+^  eV. 

Thus    the  wavelength    considerations    require    E     «    1.45   x    10      eV     if 

our    classical    approximation    is    to   be   valid. 

Consider    now   the   restrictions    imposed  by    equation    ( 43) •      We 

require   that:  .2 

D  »   


8MHS2 
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i.2  2 
h    E 

or  D  » 


8Mu(2.4)2e/f 

H 


+   .       , 

The  binding  energy  of  H   is  2.6  eV .   Thus  we  get 


£1«f%^    (*tK)  (53) 


3-93  x  10  9  Ergs 

2.46  x  10+3  eV. 

Hence    equations    ( 39)    and    ( 43)    ^re   satisfied    if   the   energy    of 
the    incident   proton   beam    is    less    than    1,000    electron    volts.       Thus 
the   classical    description    of  process    (2)    is    valid    in    the   75    to 
200    eV   range  where   the   experiment   was   performed. 
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III.   EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  used  in  this  experiment  is  the 
result  of  the  work  of  many  people.   The  overall  apparatus,  as 
shown  in  Fig.  5S  was  assembled  by  Bush  [5]>  who  demonstrated  its 
satisfactory  operation. 

A  reasonably  mono-energetic  beam  of  hydrogen  ions  ( AE  ~  2  eV ) 
is  produced  in  the  Duoplasmatr on ,  which  was  built  and  investigated 
by  Carter  [6].   The  beam  is  mass  analyzed  in  the  mass  spectrometer 
constructed  by  Strohshal  [7],  and  H   ions  proceed  on  to  the  scat- 
tering chamber  which  contains  the  target  gas,  CH,  .   The  scattering 

cell  is  placed  on  the  magnetic  axis  of  a  large  cylindrical  magnet 

+ 
at  the  position  of  maximum  field,   H_  ions,  produced  by  the  col- 

+ 
lision  of  H   with  CH, ,  exit  from  the  scattering  cell  and  are 

focused  by  the  axially  symmetric  non-uniform  field  of  the  focusing 

magnet  to  a  point  on  the  magnetic  field  axis,  a  distance  Z   from 

the  scattering  cell.   The  detector  which  is  located  at  Z  ,  is  an 

o 
/: 

electron   multiplier    capable   of   gains    up    to    10    . 

A.   THE  DUOPLASMATRON 

The   duoplasmatr on    is    shown    in    Fig.    6.       With    hydrogen    in    the 
source   and  filament   power    at    35  watts,    one   has    about     .1   amp    of 
electron    emission    from    the    filament.      An   arc    voltage   of    400    volts 
between    the    filament   and    the   Z-electrode   accelerates    these    elec- 
trons   through    the   hydrogen    gas,    ionizing    the   hydrogen    and   forming 
a   plasma.       Once    the  plasma    has    been    formed,    an    arc    current    of 
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1    ampere    can    be   sustained  between    the   filament    and    the   Z-electrode 
by   a    voltage    of   about    250    volts. 

To   extract    the  plasma,     the  positive   arc    voltage    is    switched 
from    the   Z-electrode   to    the   anode,     (see   Fig.    7)    which    is    itself  at 
a    voltage   E   above   ground.       The   energy    of    the   hydrogen    ion    beam    thus 
extracted  will    be    (Ee)    for    singly    charged    ions. 

To    enhance    the   beam    intensity   a    simple   accel-decel    system    is 
employed.       Figure   7    shows    lens    1    at    some   arbitrary    negative  po- 
tential   V.      Thus    the    hydrogen    ions    are   actually    extracted   from 
the    duoplasmatr on    at    an    energy    (V+E)e.      A    "sock"    is    fitted   over 
the    focusing   and   deflector    lenses    and   attached    to    lens    1    so   as    to 
float   at    this    negative  potential   V,    and    thus    prevent    the   beam   from 
"seeing"    ground  potential    in    this    region.      Thus    the  beam    of    hydro- 
gen   ions   passes    through    the    focusing    and   deflector    lenses    at   an 
energy    (V+E)e.      Once   the   beam   leaves    the   region    of    the   sock    it 
"sees"    ground  potential    on    the  walls    of    the   vacuum   chamber    and 
slows    to   an    energy    Ee   and  passes    on    to    the  mass    analyzer . 

This    simple   accel-decel    system    increased   the   beam    intensity    by 
a    factor    of   40 .       Thus    we   were   able   to    obtain   an   analyzed   beam   of 

100    eV   protons    (with    diameter     .5   cm   at    the   scattering    cell    80   cm 

-7 
from   the    source)    which   was    slightly    greater    than    1    x    10         amps. 

Carter    [6]    has    shown    that    the    energy    spread   of    the  proton    beam    is 

less    than    2.5%. 
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B.  THE  MASS  ANALYZER 

The  hydrogen  ion  beam  is  mass  analyzed  in  the  analyzer  con- 
structed and  calibrated  by  Strohsahl  [7].   It  employs  a  45 

bending  angle  (see  Fig.  8)  and  has  a  resolution  of  —  =  40  with 

Am 

beam  transmission  of  about  90%= 

Figure  9  shows  the  current  through  the  mass  analyzer  magnet 
that  will  pass  any  particular  species  in  the  energy  region  from 

50  to  650  volts.   Hence  .45  amps  allows  a  100  eV  beam  of  protons 

+       + 
to  travel  into  the  scattering  cell,  whereas  the  H„  and  H   are  not 

bent  enough  by  the  magnetic  field  to  be  transmitted  (see  Fig.  5). 

The  mass  analyzer  is  carefully  aligned  with  the  scattering 

cell  and  the  detector,  all  of  which  are  positioned  on  the  axis 

of  the  focusing  magnet.   (See  section  F  below.) 

C.  THE  SCATTERING  CELL 

The  scattering  cell,  shown  in  Fig.  10,  was  designed  and 
tested  by  Bush  [l].   The  target  gas,  methane,   is  bled  into  the 
scattering  cell  and  is  accurately  controlled  by  a  Variable  Leak 

valve  constructed  by  Granville  Philips.   The  pressure  of  the  CH, 

-3 
gas  in  the  scattering  chamber  (approximately  10  '  torr)  is  mea- 
sured by  the  VG  1A  ion  gauge.   The  gas  is  of  research  grade, 
(purity  99-65%) . 

The  incident  proton  beam  passes  first  through  an  adjustable 

aperature  (generally  set  at  . 5  cm  diameter)  and  then  into  the 

+  + 

scattering  chamber.   Here  the  H   reacts  with  the  CH,  to  form  H„ 

4  2 

whose  escape  from  the  scattering  cell  is  limited  by  the  geometry 
of  the  exit  slit.   The  exit  slit  is  an  annular  opening  of  36O  , 
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reduced  by  support  legs  to  about  290  ,  around  the  scattering  cell, 
separating  the  front  from  the  rear  of  the  scattering  cell.   The 
exit  wall  of  the  front  section  of  the  scattering  cell  slopes  up- 
wards at  49   and  the  exit  wall  of  the  rear  section  of  the  scat- 
tering cell  slopes  upward  at  36  .   Appendix  III  of  Ref.  5  shows 
how  these  two  angles  determine  the  target  thickness.   By  varying 
the  separation  of  the  front  and  rear  sections  of  the  scattering 
cell,  the  target  thickness  can  be  increased  from  0  to  0.8  cm. 

The  ions  collected  by  the  front  and  rear  sections  of  the 
scattering  cell  are  measured  by  separate  ammeters.   The  beam  col- 
lector collects  the  unscattered  beam  that  passes  through  the 
scattering  chamber.   This  current,  plus  the  current  from  the  rear 
section  of  the  scattering  cell,  is  recorded  as  the  total  proton 
current  incident  on  the  CH,  target  molecules. 

Three  grids  on  the  beam  collector  are  given  potentials  to 
suppress  secondary  electron  emission  and  prevents  slow  background 
ions  from  reaching  the  collector  plate.   The  beam  collector  and 
its  grids  are  placed  inside  a  metal  cup  which  shields  them  from 
other  charged  particles  and  also  acts  as  a  gas  cap  over  the  rear 
of  the  scattering  chamber. 

The  beam  collector  assembly  is  mounted  on  a  rod  which  passes  to 
the  outside  of  the  vacuum  system,  allowing  the  assembly  to  be 
pulled  away  from  the  rear  section  of  the  scattering  cell  and  swung 
out  of  the  path  of  the  incident  ion  beam.   This  allows  alignment 
of  the  entire  apparatus  and  also  allows  the  incident  beam  to  fall 
directly  on  the  detector  when  desired  for  calibration  purposes. 
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D.   THE  FOCUSING  MAGNET 

The  magnetic  field  of  the  focus ing  magnet  is  nonuniform  and 
axially  symmetric  and  has  the  focusing  property  of  bending  the  tra 
jectories  of  charged  particles  so  they  follow  a  path  similar  to 
that  shown  in  Fig .  11. 

If  one  considers  the  simpler  case  of  an  axially  symmetric 
uniform  magnetic  field,  it  can  be  shown  that  the  trajectory  of  a 
charged  particle  is  a  helix  (see  Ref.  5)»   The  particle  will  cross 
the  field  line  on  which  it  originated  after  one  complete  revolu- 
tion.  Then,  as  shown  in  Fig.  11,  particles  of  charge  q  leaving 
the  source  with  momentum  p  =  mv  at  an  angle  0  will  cross  the 
magnetic  axis  at  a  distance  Z   given  by; 

2n  p  c  os  0  ■  . 

Zo=— % t1' 

where  B  is  the  magnetic  field  strength.   Thus  a  detector  located 

at  Z  will  only  sense  these  ions  scattered  out  of  the  scattering 
o 

cell  placed  at  Z  =  0  which  have  momentum  p,  charge  q,  and  scat- 
tering angle  Q.   If  one  now  holds  the  detector  fixed  at  Z   and 

o 

increases  B  then  particles  of  larger  p  cos  0  values  will  strike 
the  detector.   Alternately  by  holding  B  fixed  and  varying  Z  we 
can  detect  particles  scattered  from  the  source  at  different  values 
of  p  cos  0.   Since  equation  1  is  independent  of  azimuthal  angle  0, 

the  entire  36O   can  be  observed  at  Z  .   For  most  conventional  scat- 

o 
ter ing  experiments  iO  •«  1  ,   Hence,  the  solid  angle  is  enhanced  by 

a  factor  of  approximately  300° 

In  our  experimental  apparatus  the  magnetic  field  is  NOT  uni- 
form, as  it  is  produced  by  a  "thin"  solenoid  shown  in  Fig.  5. 
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This  complicates  the  trajectory  equations  and  makes  numerical 

inter grat ions  necessary.   The  results  of  numerical  inter gr at  ions 

+ 
for  the  H„  ion  are  shown  in  Fig.  12.   However,  equation  1  the 

result  obtained  for  the  uniform  field  case  is  qualitatively  the 

same  as  in  our  NON-uniform  field  configuration. 

The  magnetic  field  of  the  focusing  magnet  was  measured  by 

Kelly  [8].   Figure  13  shows  the  axial  component  of  the  magnetic 

field  on  the  magnetic  axis.   Figure  14  shows  the  values  of  the 

axial  component  of  the  magnetic  field  measured  at  various  r  and 

0   values  at  the  center  of  the  coil.   Figure  15a  and  b  show  the 

radial  components  of  the  magnetic  field  as  measured  with  a  Hall 

Probe(a"t  Z   =34  cm)  across  the  geometric  axis  of  the  coil,  both 

vertically  and  horizontally.   This  demonstrates  that  the  magnetic 

and  geometric  axis  of  the  magnetic  coil  are  aligned  to  within 

0.1%.   This  allows  use  of  the  coil  geometric  axis  and  the  coil 

center -line  as  reference  for  alignment  of  the  entire  system  along 

the  magnetic  field  axis. 

E.   THE  SCATTERED  ION  DETECTOR 

The  beam  detector  and  amplifier  is  a  Model  306  Magnetic  Elec- 
tron Multiplier  whose  schematic  is  shown  in  Fig.  l6 .   The  multi- 
plier is  connected  to  a  Bendix  1122  Power  Supply  as  shown  in 
Fig.  17.   The  output  of  the  multiplier  is  fed  to  the  input  head 
of  a  Keithly  640  Vibrating  Capaciter  Electrometer . 

The  Multiplier,  developed  by  Goodrich  and  Wiley,  provides  re- 
producible current  gains  up  to  10   when  operated  at  the  potentials 
shown  on  Fig.  17.   When  operating  as  the  detector  of  scattered 
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ions  in  our  experimental  apparatus,  the  multiplier  has  to  be  oper- 
ated in  a  magnetic  field  of  up  to  200  gauss.   The  effect  of  this 
magnetic  field  is  to  reduce  the  gain  of  the  multiplier;  increasing 
field,  causes  the  gain  to  decrease. 

The  multiplier  is  operated  inside  an  aluminum  shield  to  pro- 
tect it  from  random  background.   The  various  strip  potentials  were 
chosen  to  give  the  highest  gain  and  the  lowest  noise.   The  multi- 
plier gain  was  measured  at  different  magnetic  field  values  in  the 
manner  listed  in  Appendix  I  and  fitted  by  a  12th  order  polynomial. 
See  Fig.  l8 .   The  fitted  polynomial  was  used  in  the  actual  cross 
section  calculations. 

Bush,  has  shown  that  the  gain  does  not  vary  with  gas  pressure 

-4 
up  to  2.0  x  10    torr ,  a  pressure  much  higher  than  the  multiplier 

ever  encounters  when  data  is  being  taken.   The  gain  is  independant 

of  the  incident  ion  mass  and  of  the  incident  current,  provided 

the  anode  current  is  below  its  saturation  value  of  10    amperes. 

To  prevent  slow  ions  from  reaching  the  detector,  a  grid  is 

placed  over  the  detector  shield  aperature  as  shown  in  Fig.  19. 

The  grid  also  permits  a  crude  analysis  of  the  scattered  beam  energy. 

See  section  A,  Chapter  IV. 

F.   SYSTEM  ALIGNMENT 

In  order  to  detect  particles  scattered  at  a  particular  angle 
the  detector  must  be  placed  on  the  focusing  magnetic  field  axis  at 
the  point  Z  where  the  scattered  particles  focus  on  the  axis.   The 
program  SOLANG  which  calculates  the  intercept  distance  assumes 
that  the  scattering  took  place  on  the  axis  at  the  center  of  the 
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focusing  magnetic  field.   Therefore  accuracy  in  the  angle  of  scat- 
ter determination  depends  on  how  well  the  scattering  center  and 
the  detector  are  aligned  to  the  magnetic  field. 

Figure  15  shows  that  the  geometric  center  line  of  the  focusing 
magnetic  coil  is  collinear  with  the  magnetic  field  axis  to  within 
0.1%.   This  geometric  center  line  was  used  as  the  reference  for 
aligning  all  components  to  the  magnetic  field  axis.   To  perform 
the  alignment,  cross  hairs  were  placed  on  either  side  of  the  mag- 
netic spool  to  mark  its  axis.   Then  a  small  laboratory  laser  was 
positioned  so  its  beam  lay  on  the  magnetic  axis  established  by  the 
two  crosshairs.   The  crosshairs  were  then  removed  and  the  line 
established  by  the  laser  beam  was  used  as  the  reference  for 
aligning  all  other  components  to  the  focusing  magnet  axis.   This 
method  made  possible  the  alignment  of  the  scattering  cell,  the 
detector,  and  the  focusing  lenses  to  the  magnetic  field  axis.   The 
mass  analyzer  was  constructed  by  Strohsahl,  (see  Ref.  7)  so  that 
the  ion  beam,  after  analysis,  travels  down  the  geometric  axis. 
Hence  attaching  crosshairs  to  the  front  and  rear  of  the  mass 
analyzer,  it  was  adjusted  into  alignment  with  the  laser  beam. 
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IV.   THE  EQUATIONS  DESCRIBING  THE  MOTION  OF  THE 
H*  ION  AND  THEIR  SOLUTION 


A.   THE  TRAJECTORY  EQUATION 

+ 
The  equation  of  motion  of  the  H„  ion  in  the  nonuniform  axially 

symmetric  magnetic  field  of  the  focusing  magnet,  (equation  (5)), 
is  expressed  in  terms  of  the  magnetic  vector  potential  A.   There- 
fore the  measured  magnetic  field  B  of  the  focusing  magnet  must  be 
related  to  this  quantity  before  the  trajectory  equation  can  be 
solved.   The  expansion  of  A  in  terms  of  the  field  B  is  shown  by 
Bush  [ 51  to  be: 

.  , .n      2n+l  9   B  (r=0,z) 
V^>  =LonTblTT(2)      "— ^ (D 

Ar(r,z)  =  0  (2) 

Az(r,z)  =  0  (3) 

Gagliano   [9]    fitted   the   actual    magnetic    field  measurements    along 
the   r    =    0   axis    with   a    12th    order    polynomial    in    z    such    that: 

B2(r=0,z)     =   ±     £j      vn-l  (4) 

where  I  is  the  current  through  the  magnetic  coil.   Hence  we  know 
B  (r  =  z)  and  Art(r,z)  for  any  given  focusing  magnet  current  I. 
Gagliano,  and  subsequently,  Bush,  has  shown  that  the  equation  of 
motion  of  a  particle  of  charge  q  and  mass  m  with  momentum  p  in  our 
field  configuration  is  given  by: 
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(6) 


To   solve   this    second  order    nonlinear    differential    equation   by 
means    of    the   digital    computer,    it   must    be   reexpressed  as    two   first 
order    differential    equations.      This    is    done   by   rewriting   equation 
(5)    as    follows: 


d_    /dr_ 
dz       dz 


dz  [A0  &z    j        Vdz/   lA0  dz 


l0 


a  a, 
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,2    a2 
k   "A0 


v0 
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2  ""2 
k      '  A0 
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Define: 
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V0 


Then    equation    (7)    becomes 
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dz 


[DERY(l)]     =    DERY(2) 
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$j-    [DERY(2)]     =    DERY(2)    x    FUNC(4)    +     [dERY(2)]3x    FUNC(4)       (9 
-    FUNC(5)     -     [DERY(2)]2x    FUNC(5) 

We  now  have  two  first  order  nonlinear  differential  equations 
which   can    be    solved  by    the   use    of    the    digital    computer . 

From  equation  (6)  we  see  we  need  to  know  the  momentum  of  the 
scattered  particle.  But  p  = J  2mE  where  E  is  the  energy  of  the 
scattered  particle.      Hence 


,2         2m     " f 
k      =   —   E' 

q     s 


(10) 


where  E'  is  the  energy  of  the  scattered  particle  expressed  in 
electron  volts.   Thus  if  E   is  known,  equations  (8)  and  (9)  can 
be  solved  to  yield  r  as  a  function  of  Z. 


B.   THE  ENERGY  OF  THE  H, 


To  determine  E  ,  the  energy  of  the  scattered  H_  particle, 

s  ^ 

consider  the  following  sketch  of  the  kinematics  of  the  reaction: 

H+  +  CH,  -»  H_  +  CH„. 
—       4—2      3 
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^    [DERY(2)]     =    DERY(2)    x    FUNC(4)    +     [dERY(2)]3x    FUNC(4)       (9) 
-    FUNC(5)     -     [DERY(2)j2x    FUNC(5) 

We  now  have  two  first  order  nonlinear  differential  equations 
which   can    be    solved  by    the   use    of    the    digital    computer . 

From  equation  (6)  we  see  we  need  to  know  the  momentum  of  the 
scattered  particle.  But  p  = J  2mE  where  E  is  the  energy  of  the 
scattered  particle.       Hence 

k2    J   22  E«  (10) 

q        s  v        ' 

where  E'  is  the  energy  of  the  scattered  particle  expressed  in 
electron  volts.   Thus  if  E   is  known,  equations  (8)  and  (9)  can 
be  solved  to  yield  r  as  a  function  of  Z. 

B.   THE  ENERGY  OF  THE  H+ 
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To  determine  E  ,  the  energy  of  the  scattered  H„  particle, 

s  .  2  r 
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represents    the  mass    of   the  CH„,    the  remaining  molecule.       Ignoring 

the  binding    energies,    we   can   write   the    following: 

From    the    conservation    of    energy    we   have: 

2  2  2 

%  m   v      =^mv      +  h  m   v  (11) 

oossrr  v/ 

From   conservation    of  momentum   we   have: 

m    v      =   m    v      cos    9    +   m    v      cos    0  (12) 

o   o  s    s  r    r  v       ' 

0    =   m   v      sin   0    -   m   v      sin    0  (13) 

s    s  r    r  v       ' 

We   wish   to   solve   these   three   equations    to  yield   the  relation: 

E      =    E(E,m,m,m,9) 
s  sv    o        r         s         o 

Squaring  (12)  and  (13)  and  adding,  we  eliminate  0  and  get: 

2    2  2    2  2    2 

m   v      -    2m   v   m   v      cos    9    +  m   v      =   m   v 
oo  ooss  ss  rr 

2    2  2  2 

But   from    (11)      mv      =mmv      -mmv 
v/        rr  roo  rss 

Thus    we   have   on   substituting    this    into   the  above   equation: 

2    2  2    2  2  2 

m   v       -    2m   v    m   v      cos    9+mv       -mmv      +mmv      =0 
oo  ooss  ss  roo  rss 


or  2m   E      -    4  \|  m   E  m    E    '    cos    9    +    2mE      -2mE      +2mE      =0 

—  ooyooss  ss  ro  rs 

Thus    we  get    on    rearranging: 

E(m      +m)+E(m      -m)=2   \|m   m    E   E   '     cos    9  (l4) 

svs  r '  ovo  r '  ^osos  v/ 

Now    let:       a    =   m      +   m      and   b    =   m      -   m  (15) 

s  r  or 

Squaring    (l4)    we   get; 

a2E2    +    f2ab   E      -    4   m   m    E      cos    0      E      +   b      E      =0 
sL  o  oso  Js  o 
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Solving   for    E      we    obtain: 
s 


2  2      2  2    2      2 

4   m   m    E    cos    9    -    2abE      +        V  (2abE      -    4   m   m    E      cos    0)     -    4a    b   E 
oso  o—        \|v  o  oso  '  o 


2a 


As    the   +   root    only    yields   physical    results,    in    our    case   we   have, 


on    substituting    equation    (15)    for    a    and   b: 


E      =    E 
s  o 


2mm      cos9+(m-m)(m+m) 
o   s r       o'  v    s         r  ' 

(m      +   m    ) 
v    s  r  ' 


o    o  7  o 

>  4   m  m      cos    0    -    4(m    +m    ) ( m    -m    )m   m      cos    0 
N         os  vsr/vor/os 


( m      +   m    ) 
x    s  x  ' 


Now :         m      =   m 

r  + 

CH3 


m      =    m 


Hn 


m      =   m         -    2m 
s  „+  o 

H2 


Thus    the   above    equation    simplifies    to: 


E      =   E 
s  o 


2  2  2  2 

m      +   m  m      -    2m      +   4m      cos    0 
r  r    o  o  o 


(m      +    2m    ) 
v    r  o' 


(16 


2    xi^  m      cos    9  \  m      +mm      -2m      +2m      cos    9 
*J         o  X    x r    o o o 


(m       +    2m    ) 
v    r  o' 


Thus    selecting       0  and  E      we   can   calculate   E    ,    the    energy    of    the 

scattered   H    ,    from    equation    (16). 

E 
s 
Figure   20,    a   plot    of    the   ratio  — —      versus    the    scattering    angle 

o 

0,  shows  that  in  the  10   range  from  4l°  to  51  ,  E   changes  by  less 

than  3%. 
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C.   THE  METHOD  OF  SOLUTION 

Assume  protons  of  energy  E   and  a  focusing  magnetic  current  I 

o  s»    »  0 

to  have  been  selected.   Then  equation  (4)  yields  the  magnetic 

field  B  (r  =  0,z).   Equation  ( 16 )  yields  the  energy  of  the  H   as 

2 
a  function  of  the  scattering  angle,  Q.   From  equation  (10)  k   is 

then  known  for  the  H   and  we  are  ready  to  solve  equations  (8)  and 

(9).   The  solution  is  a  trajectory  similar  to  the  ones  shown  in 

Fig.  7.   The  distance  Z  where  r  is  zero  is  recorded.   This  is 

+ 
where  the  detector  must  be  located  to  detect  H„  particles  scat- 
tered out  of  the  scattering  cell  at  the  angle  0  when  the  focusing 
magnet  current  is  set  at  I   and  the  energy  of  the  incident  beam 

of  protons  is  E  . 

r  o 

If  the  detector  distance  Z   is  changed  and  I   and  E   are  held 

o  o       o 

+ 
constant  then  H„  particles  scattered  out  of  the  scattering  cell 

at  different  scattering  angles  9  are  detected. 

Alternately,  one  can  hold  E   and  Z   constant  and  vary  I  .  This 

o       o  o 

also  allows  particles  of  different  scattering  angles  to  be  detected, 


D.   THE  COMPUTER  SOLUTION 

The  computer  solution  of  equations  (8)  and  (9)  was  originally 
worked  out  by  Gagliano,  but  has  since  been  rewritten.   Program 
SOLANG,  explained  in  Appendix  II,  employs  the  Hamming  Modified 
Pr edictor -Corr ector  method  to  integrate  equations  (8)  and  (9). 
This  is  the  DHPGC  subroutine  of  the  IBM  System  36O-67  Scientific 
Subroutine  Package. 
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The  trajectory  equations  are  solved  for  a  range  of  scattering 

angles  for  a  particular  E   and  I  .   The  results  are  stored  and 

o  o 

later    used    to    compute    the   solid   angle   subtended   by    the   detector    as 
it    intercepts    each    of    these   trajectories. 

E.       THE    SOLID  ANGLE 

The   solid   angle   subtended   by   the   detector    is    given   by    the 

relation: 

dQ  (0,0)    =    sin   0   d©    d0 

But    since    our    system    is    axially    symmetric    this    becomes 

dQ  (0)       =    sin   0   d0cfd0 

290 
=    2tt  (3S0)    ^n    0    d© 

=    1  .612  tt    sin    0   d©. 

This    is    because   only   290      of    the   available   360      of    the  axial    angle 

0  are   accessible    to    the   detector    due    to    the  construction    of    the 

scattering   cell.       (See  Chapter    II,    Section   3.) 

From    the    trajectory   plots    on    Fig.    15    it    can    be   seen    that   all 

particles    scattered    into   0  will    cross    the   axis    at   Z      within    AZ    • 
r  00 

This    AZ      determines    the   angular    resolution    A©   which    the    detector 
o 

sees    and    so   allows    determination    of    the    solid   angle.       But    AZ      is 

"  o 

in    turn    determined   by    the    detector    geometry   and    the  particular 

trajectory    under    observation.      The   relation   between    A©   and    AZ      can 

J  '  o 

be    found    from    the    functional    relationship: 

Zq   =    f(E,    B,    ©,  p  ,    T],    C) 
where    E    =    energy    of    incident    ion. 
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B  =  magnetic  field 

K   =  target  width  from  system  axis  along  a  radial  line 
n  =  target  thickness  from  system  center  along  the  Z  axis 
C  =  constant  for  each  reaction  including  mass  ratios, 
inelastic  energy  losses,  etc. 
Differentiating  this  gives: 

9z  dz  dz  dz  dz  oz 

4Zo  -  sf  iE  *  si2  iB  +  do2  ie  +  df  7° +  pr  a11  +  sr 4C 

But    AB   =    0   and    AC   =   0    for    any    specific    reaction   and   measurement 

Thus,    solving    for    A©   we   get: 

fbZ  \  /dz  \  (dz 


°      VdE^    Q,e,T1  ^Wb.Q.T)     I       l*V    E>Q,p 

'dZ    V  ' 

_o 


No  analytical  expression  exists  for  the  various  differentials. 
However,  computer  trajectory  results  obtained  by  varying  E,0,v 
and  7]  show  the  effect  of  the  differentials  is  to  change  A0  by  no 
more  than  by  10%  (see  Ref.  10).   Thus  the  angular  spread  can  be 
expressed  as: 

49  "(H;)  4Zo  <17> 

Hence    the   solid   angle   becomes 

dQ  (0)    =    1.6l2TTsin   e/l^)AZ0  (l8) 

Now,    since   AZ      depends    on    the   angle   0       (see   Fig.    21   and 
o  o 

Ref.    5    for    further    details)    and  0       in    turn    depends    on    the   scat- 
tering   angle   0   we    see   that    dD  must    be   evaluated   at    each   angle   of 
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Figure  21.   Angular  and  Energy  Acceptance  of  Dector 
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scatter.   This  is  implemented  in  the  latter  subroutine  of  the 
computer  program  SOLANG . 

Computer  results  indicate  that  the  angular  resolution  given 
by  equation  (17)  is  approximately  1  . 

F.   THE  CROSS  SECTION 

Consider  a  beam  of  monoener get ic  ions  I.  incident  on  a  target 
of  N  molecules.  A  detector  located  at  an  angle  0  subtends  solid 
angle  dQ  (0)  . 


I. 

l 


dQ(0) 


«i 


N. 


The  number  of  particles  scattered  into  the  detector,  N  ,  is 

s 

usually  measured  as  a  current  I   where: 

s 

I      =    N      dQ  (0) 
s  s  v     ' 


The   cross    section    is    then    defined   as 

I 

dCT(Q)     =    - 


it  v     ' 


(19 


But  I   is  magnified  G  times  by  the  multiplier  gain.   Hence  I  , 
the  actual  detector  current  measured  is: 

I~  =  I   G 

D     s 
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Hence    (19)    becomes: 

I 
da  (0)    = 


I  .    G   N      dfi(O) 
1  t  v     ' 


Bu 


t    N    ,     the  number    of   target   particles,     is    given    by 


N      =    P       P  t 
t  sc 


Where   P         =    pressure    of    target    gas    in    scattering   cell    in    torr 


i  ft  ^ 

r  =    3.536    x    10        particles/cm      torr 


P 
t  =  target  length  in  cm 

Hence  we   have: 

dCT'e»=  1.  cp%t  dn(o)   CI"2  <20> 

1  SC  *  v       ' 

Because    of    the   construction    of    the   scattering   cell    (see 

Section    3    of   Chapter    II)     the    target    thickness    depends    on    the 

scattering   angle  0.       Bush,    has    shown    this    relation    to   be: 

2.54    sin    (36°-8)    if  o 

sin    0 

t    =    S  if   36°  <    0  <    49°       (21) 

t  =  s  -  2-°  sin   (Q-49°       if  49°  <  0 
s  in   0 

where    S    is    the   minimum   separation    of    the   front    of    the   scattering 

cell    from    the   rear    of    the   scattering   cell,     (which    is    adjustable 

from    outside    the   vacuum   system). 

Computer    Program   CSVSZ,    listed    in   Appendix    I,    calculates    the 

differential    cross    section    in    the  Lab.    and   Center    of  Mass    coordin^ 

ate   systems    using    equations    (20)    and    (21)    and    the   solid   angle 

determined   by    the  program   SOLANG . 
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V.   EXPERIMENTAL  RESULTS 

A.   ELIMINATION  OF  BACKGROUND 

+ 
When  H   collides  with  CH,  three  general  types  of  collision  may 

occur  which  can  produce  large  detector  signals. 

Event  (a)  The  rearrangement  collision  where  the  H   captures  a 

+ 
proton  to  form  H„ . 

Event  (b)  A  scattering  collision  of  the  H   with  a  hydrogen  atom. 
Event  (c)  A  scattering  collision  of  the  H   with  a  carbon  atom. 
Events  (b)  and  (c)  are  expected  to  have  cross  sections  that  are 
many  orders  of  magnitude  larger  than  those  of  event  (a).   Hence 
appropriate  steps  were  taken  to  insure  events  (b)  and  (c)  were 
not  measured  by  our  detector.   To  eliminate  events  (b)  and  (c) 
we  look  first  at  the  energy  and  momentum  of  the  scattered  parti- 
cles in  all  three  events. 

Consider  first  event  (a).   The  capture  process  may  be  approxi- 
mated by  the  reaction  shown  below: 

H+    +    CH,    ->    H+    +    CH0 

—  k        —  3 


o 


m 


v,°     Jf\  s 


S-/L    P 


o  o  m 


\ 


The    energy    of    the    incident    H     particle    is 

2 
2        P- 


E      =    %  m    v      =    , 
o  o    o         2m 


(1) 
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Now,  from  conservation  of  momentum  we  see  the  momentum  of  the 

+ 
scattered  H  particle,  represented  by  p  , related  to  p   by: 


/  e°     s 
p   ■«  p   COS  u5   =  

°     S  {2 


(2) 


where  we  assume  the  mass  of  the  CH  ,  represented  by  m  ,  to  be 

3  t 

infinitely  large.  Hence  in  the  capture  process,  the  scattered 
particle  has  momentum  p  =  ^~2~  p  .  The  energy  of  the  scattered 
particle  is  then: 


E  = 


P     P 
rs     r  o 


s    2m      m 
s      s 


But  for  the  H^  particle,  m   =  2  m 
2  r  s       o 

2 


Hence: 


o =  E 

2m      o 


(3) 


Consider    now    event    (b),    the   scattering    of   the   H      from   a 
hydrogen    atom: 


H+    +    H  -    H+    +    H 


Since   in      =    m      and  p      = 

t  s  c  t 


p      we   see   that: 
s 


Hence: 


Then: 


Po    =    2   Ps    °OS    /f5°   =  \P"     P< 
Po 


E      =    

s         2m 


4m 


=   h 


(4) 
(5) 
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+ 
Thus  the  elastically  scattered  H   from  H  has  a  momentum  that  is 

Po 
approximately   and  an  energy  that  is  approximately  \   E  • 

(F  ° 

Finally  consider  event  (c),  the  scattering  of  the  H   from  a 

carbon  atom: 

H   +  C  -»  H   +  C 


Assuming  the  carbon  to  be  infinitely  massive  we  see  that: 
Hence: 


But: 


Thus:  p^  =  p^  (7) 

+ 
Hence  The  elastically  scattered  H   from  C  has  momentum  p   and 

^o 

energy  E  . 
^    o 

Summarizing  the  results  of  this  very  approximate  analysis  in 

Table  1  we  get  for  an  incident  beam  of  protons  of  energy  E   and 

momentum  p   the  information  shown.   Note  that  the  scattered  H 
*o 

produced  in  events  (b)  and  (c)  has  momentum  less  than  the  H„ 


E        = 

s 

E 
o 

2 
Ps 

= 

2 

Po 

2m 
s 

2m 
o 

m 
s 

= 

m 
o 

Ps 

= 

Po 

produced  in  event  (a) 
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Possible   Events 

Scattered 
Particle 

Momentum    of 
Scattered  Particle 

Energy    of 
Scattered 

Par  tide 

Event    (a)    Capture 

H+    +   CH,^    H*    +    CH0 
—              4—2             3 

+ 
H2 

^   Po 

E 
o 

Event    (b)    Scatter 
H+    +    H  -    H+    +    H 

+ 
H 

$TPo 

*   Eo 

Event    (c)    Scatter 
H +    +    C  ■*    H+    +    C 

+ 
H 

Po 

E 
o 

TABLE   1.       The  Three   Possible  Scattering   Processess. 

Plots    of    the   scattered    ion    trajectories    resulting    from   each 
of    the    three    events    tabulated   are   shown    in    Fig.    22,    for    the  parti- 
cular   case    of    E      =    100    eV   and    I    =    6.0   amps.       Note    that    the    two 
o  c 

+ 
elastically    scattered  H      ions    cross    the   axis    at   a    Z   less    than    that 

+    .  + 

of    the   H      ion.      To   prevent    these   unwanted  H     particles    from 

reaching    the   detector,    a    baffle,    5   cm    in    radius    was    positioned 

over    the    beam    collector    as    shown    in    Fig.    22.      This    allowed   only 

+ 
the   H      produced   by    the   rearrangement    collision    to    reach    the 

detector . 

Figure   23   contains    plots    of    the   detector    current    versus    the 

detector    grid    voltage    for    the   case    of   E      =    100    eV   and    I    =    6.0   amps 

and  with    the    detector    located   at    Z   =    46    cm    (this    corresponds    to   a 

scattering   angle   0    of   46.92    ).       Figure   23    -(i)    was    recorded   before 

the   baffle  was    placed    in    the  path    of   the   H      ions    and   Fig.    23    -    (ii) 

was    recorded   after    the   baffle  was   positioned,    as    shown    in    Fig.  22. 
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om 


ween 


Consider  first,  Fig.  23  -  (i).   A  grid  voltage  of  25  volts 
eliminates  the  slow  ions.   As  the  detector  grid  voltage  is  varied 
from  25  to  50  volts,  we  continue  to  collect  the  ions  scattered  fr 
the  scattering  cell  due  to  processess  (a)  and  (b)  and  (c).   Bet 
60  and  90  volts  we  collect  only  these  particles  with  energies 
greater  than  90  volts;  ie,  the  ions  which  result  from  processes 
(a)  and  (c)  only. 

Figure  23-(ii)  now  shows  that  the  baffle  does  indeed  eliminate 
particles  resulting  from  processes  (b)  and  (c).   Thus  with  the 
grid  voltage  set  anywhere  between  30  and  90  volts,  and  the  baffel 
in  position,  we  only  detect  the  H„  that  results  from  the  re- 
arrangement collision. 

B.   DEMONSTRATION  OF  VALIDITY  OF  DETECTOR  SIGNAL 
From  equation  (20)  of  Chapter  IV  we  have: 

dCT<Q>=  1.  p  lDetGdQ(s)  c™2  <8> 

1         SC    x  v       ' 

If    the   angle    of   the   scatter     0  is    held   constant,    then    this    re- 
duces   to: 

I 
dcr  =   K  J-!—  (9) 

i    sc 

where  K.  =  constant. 

Consider  now  the  case  where  I.  is  held  constant.   Then  (9) 

1  v  ' 


reduces  to 


1,=  C.  P  (10) 

D     1   sc  v   ' 


wh 


ere  C,  =  constant.   Experimental  data  plotted  in  Fig.  24  obeys 


this  linear  relationship.   The  fact  that  I   does  not  extrapolate  to 
zero  is  due  to  the  background. 
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Figure  23.     Detector  current  vs.  Detector  grid  voltage 
( i  )   Without  baffel 
(ii)    With  bafftl 


SO  40 

Otttctor  grid  voltagt 
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Consider  next  the  case  where  P    is  held  constant.   Then 

sc 


equation  (9)  reduces  to; 


In  =  C_  I. 
D     2   1 


(11) 


where  I.  =  constant.   Experimental  data  plotted  in  Fig.  25  follows 
this  linear  relationship.   Thus  our  detector  signal  has  the  proper 
dependance  on  the  various  parameters  and  behaves  like  a  scattered 
particle  current. 

For  the  differential  scattering  measurements  made  with  100  eV 

incident  protons,  the  proton  current  at  the  beam  collector  was 

-7 
approximately  10    amps  and  the  current  to  the  rear  of  the  scat- 

_o 

tering  cell  was  about  10    amps.   The  methane  gas  pressure  was 

-4  -3 

varied  from  1  x  10    torr  to  3  x  10    torr .   The  detector  signal 

-12 

under  these  operating  conditions  varied  from  1  x  10     amps  to 

3  x  10     amps.   This  represents  a  flux  of  a  few  hundred  particles 

per  second.   With  no  gas  in  the  scattering  chamber  the  detector 

-13 
current  fell  to  about  5-0  x  10     amps. 

The  UGIA  Ion  Gauge  used  to  measure  the  methane  pressure  in  the 

scattering  cell  was  calibrated  using  a  Capacitance  Manometer.   The 

relationship  between  the  pressure  measured  by  the  ion  gauge 

(P        ,)  and  the  actual  pressure  which  was  measured  by  the  mano- 
v  measured' 

meter  (P      , )  is  linear  in  the  region  of  interest  as  is  shown 
v  actual 


below 


actual    x    10 
torr 


slope    =    0.625 


Figure   24.       Detector    Current    versus    Pressure    of   Target    Gas 
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Figure  25.    Detector    Current   versus    Incident   Current. 
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C.   H   +  CH,  SCATTERING  DATA 

The  cross  section  for  the    formation  of  H   was  measured  at 
different  angles  of  scatter  by  varying  the  detector  distance  Z 
while  holding  the  magnet  current  I  fixed.   As  pointed  out  in 
Section  B  of  Chapter  II  the  cross  section  is  expected  to  be  a  pro- 
nounced peak  at  46 . 9  .   Hence  the  detector  was  swept  through  the 
range  of  scattering  angles  from  43   to  49  •   This  was  repeated  for 
various  target  gas  pressures  for  each  value  of  the  energy  of  the 
incident  proton  beam.   The  energies  investigated  were  70,  85,  100, 
150  and  200  electron  volts.   The  scattering  data  recorded  at  each 
of  these  energies  is  shown  in  Figs.  26,  28,  30 ,  32  and  33 
respectively.   A  prominent  peak  exists  in  the  70,  85,  and  100  eV 
data.   The  150  and  the  200  eV  data  shows  no  apparent  peak  in  the 
cross  section . 

The  data  was  analyzed  in  the  following  manner.   From  the  cross 
section  measurements  made  at  each  energy  an  average  value  of  o ( 0 ) 
was  computed  and  sketched  in  as  a  solid  line.   The  assumed  back- 
ground was  then  drawn  as  a  dotted  line.   Subtracting  out  the  back- 
ground from  the  70,  85,  and  100  eV  data,  the  points  connected  by 
the  solid  curve  in  Figs.  27,  29,  and  31  were  obtained.   To  find 
out  what  the  actual  peak  in  the  scattering  data  looked  like  before 
it  was  detected  by  our  detector  of  approximately  1   resolution, 
the  effect  of  the  angular  width  of  the  detector  was  determined  and 
the  actual  curve  unfolded.   The  curves  dotted  on  Figs.  27,  29,  and 
31  represent  the  unfolded  original  peak  in  the  scattering  data. 
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Figure   26.    Cross    section   for 
H+    +   CH,   -•    H*    +   CH 

at    70    ev.    (background 
included) 
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Figure   27.    Cross    section    for 


H2    +   CH4 


at  70  ev.  (background 

subtracted  out) 
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Figure   28.    Cross    section    for 

H+    t   CH,    "*    H_    +    CH_ 
k  2  3 

at  85  ev  (background 
included) 


Figure  29.  Cross  section  for 


H+  +  CH4  -  H+  +  CH3 


at  85  ev.  (background 

subtracted  out) 


Figure   30.    Cross   section   for 


H+   +  CH 


H2    +  CH3 


at  100  ev  (background 
included) 
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Figure   31.    Cross    section   for 

H+    +   CH^  -    H*    +  CH. 


at   E  =   100   ev 

(background  subtracted  out) 
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Figure   32.    Cross   section   for 

H+    +  CH^.-    H*   +  CH3 

at  150  ev  (background 
included) 
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Figure   33-    Cross   section   for 

H+    +  CH,    -    H*    +  CH, 


12.0 


at    200   ev    (background 
included) 
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These  peaks  are  then  the  actual  cross  sections  for  the  rearrange- 
ment reaction: 

H+  +  CH,  -•  H*  +  CH„ 
-      H    -2      3 

+ 
as    the   energy    of    the  H      varies    from   70    to    100    eV .       All    the   peaks 

exhibit    the   same   characteristics:    a    sharp    leading    edge  at    the    lower 

angles,    a    flat    top    and  a    long    decreasing    tail    extending   out    to 

about    51    •      The  peaks    dependence    on    energy    is    exhibited    in    the 

following    table: 
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70   eV 

-22       2 
5.95x10         cm 

45.15° 

1.7° 

10.2 

!  '" 

20.6 

45.87° 

85    eV 

-2?        2 

4.26x10         cm 

45.1° 

1.5° 

6.56 

11.56 

45.6° 

100    eV 

-22       2 
3.0x10         cm 

46.25° 

1.5° 

4.62 

7.6 

46.7° 

^150    eV 

-22       2 
^1.0    xlO         cm 

Not 

Not 

Not 

Not 

Not 

Observed 

Observed 

Observed 

Observed 

Observed 

Table  II.   Summary  of  Experimental  Data 

The  peak  does  exhibit  the  expected  physical  properties.   As  the 

energy  increases  the  magnitude  of  the  peak  decreases,  the  width 

.  ,   o 
narrows  and  the  angle  at  which  the  peak  occurs  moves  towards  46.9 
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The  energy  dependence  of  the  experimental  cross  section  is  compared 
to  the  theoretical  cross  section  in  Fig.  34.   Curve  (a)  represents 
o"(0)    ,  ie.  the  magnitude  of  the  flat  top.   Curve  (b)  represents 
the  area  under  the  flat  top.   Curve  (c)  represents  the  total  area 
under  the  curve. 
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Figure  34-  Energy  dependence  of 
cross  section 
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V.   CONCLUSION 

The  close  agreement  between  the  experimental  data  and  the 
predictions  of  the  ion  molecule  rearrangement  theory  proposed  by 
Bates,  Cook  and  Smith  is  evident  in  the  following  observations: 

(1)  The  sharp  peak  in  the  cross  section  predicted  by  Bates  et.al., 
is  observed  in  the  70,  85s  and  100  eV  data,,   The  fact  that  the 
peak  is  not  observed  on  the  150  and  200  eV  data  is  not  sur- 
prising for  extrapolation  of  the  experimental  data  to  150  eV  shows 

-22   2 
that  the  magnitude  of  the  cross  section  should  be  1.0  x  10    cm  . 

The  150  eV  experimental  data  on  Fig.  32  shows  that  a  peak  of  this 

magnitude  or  smaller  cannot  be  observed. 

(2)  The  theory  predicted  the  peak  to  occur  at  46.9  .   The  peaks 
observed  were  located  between  45°1   and  46.25  ,  the  latter  cor- 
responding to  the  higher  energy.   As  the  accuracy  of  the  approxi- 
mations made  in  the  theory  increases  with  higher  energy,  one 
expects  the  location  of  the  peak  in  the  experimental  data  to  ap- 
proach 46.9   as  the  energy  is  increased.   This  trend  is  observed 
exper  i men tally . 

(3)  As  is  evident  from  Fig.  34s  "the  experimental  cross  section  has 
an  energy  dependence  similar  to  that  predicted  by  the  theory. 

(4)  The  cross  section  predicted  by  Bates  et.al.,  is  an  upper  limit. 
Hence  it  is  expected  that  the  experimental  cross  section  might  be 
smaller  than  the  theoretical  cross  section.   The  experimental  cross 
section  shown  in  Figt  34  is  in  fact  about  a  factor  of  30  smaller 
than  the  theoretical  upper  limit. 
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(5)  As  some  of  the  approximations  made  in  the  theory  may  no  longer 
be  valid  at  low  energies  one  expects  the  peak  to  broaden  as  the 
energy  decreases „   This  trend  is  indeed  observed  in  the  100,  85, 
and  70  eV  data  shown  in  Figs.  31?  29,  and  27  respectively. 

It  is  felt  that  the  facts  mentioned  above  constitute  a  satis- 
factory verification  of  the  ion  molecule  rearrangement  theory  of 

+ 
Bates,  Cook  and  Smith  as  applied  to  the  formation  of  H„  in  the 

r  eaction 

H+  +  CH,  -  H*  +  CHL  . 
-       -4-2      3 
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APPENDIX  I 


MULTIPLIER  GAIN  MEASUREMENT 

The  following  sketch  shows  how,  in  principle,  the  gain  measure- 
ment can  be  made. 


Dynode 
Str  ip 


By  recording  the  incident  current  I.  and  the  detector  output  I 

the  gain  can  be  computed  from: 

I 


G  = 


D 


To  measure  I.  the  Dynode  Strip  was  used  as  a  Faraday  Cup.   Two 

problems  then  arose: 

o 

(1)  The  Dynode  Strip,  because  of  its  large  resistance  (10  Q) 

retains  a  large  residual  negative  charge.   Thus  in  using  the  Dynode 
Strip  as  a  Faraday  Cup,  I.  had  to  be  greater  than  10     amps  so 
as  to  insure  that  this  residual  negative  charge  did  not  adversely 
affect  the  measurement, 

(2)  The  Multiplier  saturates  when  I   ^  5  x  10    amps.   Thus  we  must 
have  I.  <  (5  x  10"  )G. 

To  overcome  these  two  difficulties  the  gain  measurement  was 
performed  in  the  following  sequence  of  steps. 
Step  1.   A  beam  of  about  10  '   amps  was  focused  on  the  Beam 
Collector . 
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Step    2.    The    Beam   Collector    was    swung   away    from    the   back   of   the 
Scattering   Cell.       This    allowed    the    ion   beam   to    fall 
directly    on    the   multiplier.       With    F.S.I,    and    D.S.O. 
grounded,     (see   Fig.    17)    and   a   Keithley    610   Micro-Micro- 

ammeter    connected    to    the   D.S.I,    this    incident    current 

1  -10 

was    measured    to   be    I  .     (-«    10  amps). 

Step    3-    The  Multiplier    was    now   restored    to    the   configuration    shown 
in    Fig.    17-      The    Dynode  Voltage   was    reduced    to   approxi- 
mately   1,000    volts    (this    insured    that    the   anode    current    to 

-ft 
be   measured    in    Step    4   was    less    than    10        amps    the   satu- 
ration   current    of    the  Multiplier). 

t 

Step  4.  The  Multiplier  was  turned  on  and  a  detector  current  I 

-7 
(-«  10    amps)  was  measured  at  the  anode.   Thus,  with  the 

Multiplier  settings  as  selected  in  Step  3  the  gain  was 


G   =  n- 


1 


Step  5-  The  incident  be*am  I.  was  now  reduced  by  adjusting  the 

focusing  electrodes  (see  Fig.  5)  a  significant  amount  to 

i 
cause  I   to  decrease  at  least  two  orders  of  magnitude. 

D  * 

ii  -9 

This    new    detector    current    T       (^a    10        amps)    was    then    re- 

D  it 

I 

corded.      Thus    the    incident    beam   must    have    been    I.    =   — - 1 —    , 

l         G 

Step    6.    The    Dynode  Voltage   was    turned   up    to    1,500    volts.       The 
Multiplier    now    had    the   configuration    desired  for    the 
experimental    measurements    to   be   made    later    on.       The   de- 
tector   signal    I      now   recorded   allowed    the   gain    of    the 
Multiplier    in    this    configuration    to   be    computed   from: 
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The  above  procedure  was  repeated  at  various  values  of  detector 
distance  Z  and  focusing  magnet  current  I  selected  so  as  to  give  a 
measure  of  G  versus  the  axial  magnetic  field  B.   The  results  of 
such  measurements  is  shown  in  Fig.  l8. 
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APPENDIX  II 

The  two  computer  programs  SOLANG  and  CROSEC  are  listed  in  the 
following  chapter . 

SOLANG,  is  the  program  which  integrates  the  equation   of  the 

motion  of  the  scattered  H   ion  ( equations  (o)  and  (9)  of  Chapter  IV) 

for  predetermined  values  of  the  Focusing  Magnet  Current  I,  and 

the  Incident  Proton  Energy  E  over  a  range  of  scattering  angles  0. 

The  integration  is  preformed  using  the  D.H.P.G.C.  Subroutine  of 

the  IBM  System  36O-67  Scientific  Subroutine  Package  in  the  MAIN 

portion  of  the  program,  and  in  the  two  Subroutines  FCT  and  OUTP . 

The  resulting  distances  Z  ,  where  the  trajectory  crosses  the 

o 

magnetic    axis    and   the    corresponding    angles    9      (see    Fig.    21)    are 
stored   and   later     tabulated. 

The  program   now    continues    on    to    the   Subroutine    DOMEGA   and    the 

succeeding    subroutines,    to   compute    the   solid   angle   dQ(Q)     for    each 

+ 
of    the  H      trajectories.       The   solid   angle    is    computed   using 

equation    (l8)    of  Chapter    IV,    which    is: 

d0(9)  =    1.612H  sin    0  (  ~-\  AZ 

^ZoJE,I  ° 

AZ      has    been    determined   by    Bush,     ( Ref .    5)    and    is    given    by    the 

following    relations;     (see    Fig.   (19)    for    explanation    of    symbols 


used) 


I    for       32.2°  >    9      >    26.6° 

AZ         =    (f+a)    cotan    9       -    2b 
o  v  '  o 
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II   for      26.6°  *    0     >    17.3° 

AZ      =    (h+a)    cotan   9       -    b 
o         v  '  o 

III    for       17-3°  s      Q      >    o° 

o 

AZ      =    2h    cotan   0 
o  o 

The   relation/  c— — \  is    determined   by    fitting    the  results    of   Z 

V.        o  )  E,  I 

versus    0    that   were   obtained   in    the   trajectory    integrations    by   a 

azo 

polynomial  Z   =  f(0).   We  then  evaluate  — r—   at  the  various  0 
^  o     v  '  $0 

values.   In  this  manner  the  solid  angle  dQ  is  computed  for  various 

values  of  0. 

In  order  that  the  results  may  be  tabulated  for  various  Z  values 

it  is  necessary  to  fit  the  results  of  0  versus  Z   that  were  ob- 

o 

tained  in  the  trajectory  integrations  by  a  polynomial  9  =  f(Z). 
Hence,  on  choosing  a  Z,  the  corresponding  0  can  be  computed  and 
in  the  manner  described  above,  the  solid  angle  can  be  determined. 

CROSEC  is  the  program  which  computes  the  cross  section  from 
the  experimental  measurements  of  the  Detector  Current  I   and  the 

Incident  Proton  Current  I.  which  were  recorded  at  various  detector 

l 

distances  Z  for  a  particular  proton  energy  E   methane  pressure  P 

and  a  particular  Focusing  Magnet  Current  B 

r  '  mag 

Because    the    data    is    recorded  at    various    detector    distances    Z, 

the   coefficients    of    the   polynomial    0    =    f(Z)    that    were   determined   in 

the  program   SOLANG  are   read    into   CROSEC    on    data    cards    and    so   allow 

the   0    corresponding    to   each    Z    to   be    determined.       Rather    than    read 

into    this   program    the   solid   angle   at    each  point    Z,    the   coefficients 

of    the  polynomial    0      =    f(0)    that   were    determined    in    SOLANG   are   read 

in    on    data    cards.       Hence    using:    dG(0)  =    1.612TT  sin    0/- —  )  AZ        and 

lOZ  /        o 
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the  relations  for  AZ  referred  to  in  the  discussion  of  the  program 
SOLANG,  we  compute  the  solid  angle  at  the  0  corresponding  to  the  Z 
where    the   data   was    recorded.       Hence   using    equation    (20)    of  Chapter 

IV:  I 

da(Q>  =    I.P        et    G   dQ(9) 
1    sc   v  s     ' 

we  compute  the  cross  sections. 

Using  the  transformation  from  the  Laboratory  to  the  Center  of 

+  + 

Mass  coordinates  for  the  reaction  H   +  CH,  ~*  H_  +  CH„  we  can  de- 

-      ^-23 

termine  9„   _  „  and  hence  a(0)_   _  w.   The  results  are  both 

C    of   M  *       C   of  M 

tabulated   and  plotted. 
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